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ABSTRACT 

Context. The study of abundances in the nucleus of active galaxies allows us to investigate the evolution of abundance by comparing 
local and higher redshift galaxies. However, the methods used so far have substantial drawbacks or rather large uncertainties. Some 
of the measurements are at odds with the initial mass function derived from the older stellar population of local elliptical galaxies. 
Aims. We determine accurate and reliable abundances of C, N, Ne, and Fe relative to O from the narrow absorption lines observed in 
the X-ray spectra of Mrk 509. 

Methods. We use the stacked 600 ks XMM-Newton RGS and 180 ks Chandra LETGS spectra. Thanks to simultaneous observations 
with INTEGRAL and the optical monitor on-board XMM-Newton for the RGS observations and HST-COS and Swift for the LETGS 
observations, we have an individual spectral energy distribution for each dataset. Owing to the excellent quality of the RGS spectrum, 
the ionisation structure of the absorbing gas is well constrained, allowing for a reliable abundance determination using ions over the 
whole observed range of ionisation parameters. 

Results. We find that the relative abundances are consistent with the proto-solar abundance ratios: C/O = 1. 19±0.08, N/O = 0.98±0.08, 
Ne/0 = 1.11±0.10, Mg/O = 0.68±0.16, Si/O = 1.3±0.6, Ca/O = 0.89±0.25, and Fe/O = 0.85±0.06, with the exception of S, which 
is slightly under-abundant, S/0 = 0.57±0. 14. Our results, and their implications, are discussed and compared to the results obtained 
using other techniques to derive abundances in galaxies. 

Key words, galaxies: Seyfert - quasars: individual: Mrk 509 - galaxies: active - X-rays: galaxies 



•rH 1. Introduction the plasma is difficult. In contrast, in the X-ray band a multi- 

tude of ions have observable transitions, sp anning about four 

5^ The X-ray emission from active galactic nuclei (AGN) is gen- order s of magnitude in ionisation parameter dSteenbrugge et al.1 

. , erally thought to be Comptonized emission originating from [2005I: lHolczer"itd] 12007). but the spectral resolution is signif- 
near the accretion disk surrounding the supermassive black hole, j^antly poorer and the different velocity components cannot be 
which accretes gas from the host galaxy In general, the contin- f^Hy distinguished. The most reliable results are therefore ob- 
uum IS well described by a power-law component and a soft ex- joined by performing a joint analysis of the U V and X-ray spec- 
cess. In about 50% of Seyfert 1 galaxies ( jCrenshaw et al. || 2003|) , tra, preferentially obt ained simultaneously dSteenbrugge et al.l 
narrow absorption lines from highly ionised gas are observed, |2005t lCostantini et alJl2007HCostantinill201 0^. so that the veloc- 
with the gas generally outflowmg at a velocity of between 100 jty structure obtained from the UV spectrum can be combined 
and a few 1000 km s ' (Kaast ra et al. | |2002|l . This absorption ^j^h the ionisation structure determined from the X-ray spec- 
can be studied in detail in the UV and X-ray part of the spec- trum 
trum. The spectral resolution is significantly higher in the UV, 

but only a few absorption ti-ansitions are found in this part of Although several detailed studies of warm absorbers (i.e. the 
the spectrum. Therefore, determining the ionisation structure of ionised outflow observed through narrow absorption lines in the 
X-ray spectrum) have been accomplished in the past few years. 
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the warm absorber are only now being determined. We there- 
fore executed a large observing campaign targeting the Seyfert 1 
galaxy Mrk 509 to help us resolve these questions (Kaastra et al. 
1201 lah . The absorber is generally assumed to originate either 
from the accretion disk or the obscuring torus, thus the abun- 
dances determined should closely match those of either the ac- 
cretion disk or the obscuring torus. The abundances determined 
should be representative of the abundances in the vicinity of the 
supermassive black hole and the nucleus of Mrk 509. 

In the case of Mrk 509, the soft excess originates from 
an optically thic k plasma with a temperature of 0.2 keV 
(iMehdipour et al.l l2011). while the hard X-ray power-law likely 
originates in an optically th in plasma with a much higher tem- 
perature (Palt ani et al.ll201 ll) . An absorber is detected with a v e- 
locity range between -426 and -1-219 km ' (iKriss et al.ll201ll) . 

Determining the relative abundances in the nucleus allows us 
to study the enrichment processes in the host galaxy of Mrk 509. 
Furthermore, the relative abundances that we measure can be 
compared with the (relative) abundances measured for active 
galactic nuclei (AGN) at high redshifts, to determine the abun- 
dance evolution of AGN and thus the history of the enrichment 
processes prevalent at different epochs. At high redshifts, abun- 
dances are usually determined from the relative strength of two 
or more broad emission lines shifted into the optical band. This 
broad-emissi on-line gas is emitted w ithin 0.1 pc of the central 
black hole ( Peterson & Wan del"2000^. In Mrk 509, the size of 
the broad line region (BLR) determined from H/3 measurements 
is 80 light-days (0.067 pc) (ICarone et al.l il996). The absorbing 
gas is probably located at a somewhat larger distance, although 
both aspects of the AGN phenomenon are likely to have similar 
abundances because of their close origin. Using this assumption 
allows one to compare the abundances determined by these dif- 
ferent methods, and thereby the abundances derived for different 
redshifts. 

The main advantage of using X-ray absorption lines over 
the optically detected broad emission lines is that the ionisation 
structure of the absorber can be accurately determined, which is 
crucial in differentiating the effects on the plasma of the abun- 
dances from those of the ionisation structure. Unless the ioni- 
sation structure is well known, one can only consider ions that 
have their peak ionic column density at the same ionisation pa- 
rameter to determine the abundances. There are no hydrogen 
transition lines in the X-ray spectra; only continuum absorp- 
tion is present. However, the bound-free absorption from H is 
degenerate with the continuum model used, which is not a pri- 
ori known. Therefore, only relative abundances are determined 
throughout this paper. In addition, unless mentioned otherwise, 
all the relative abundanc es are with respect to proto-solar ones 
dLodders&Palmd 12009 . see the last column in our Table [TJ. 
These abundances were determined from chondritic meteorites 
and the photospheric abundances of the Sun. Here we measure 
the abundances of C, N, Ne, Mg, Si, S, Ca, and Fe relative to 
O for gas analysed based on its narrow absorption lines in the 
stacked XMM-Newton reflection grating spectrometer (RGS) 
and Chandra low energy transmission grating (LETGS) X-ray 
spectra of Mrk 509. We thus determine the relative abundances 
of those elements that result from the different enrichment pro- 
cesses occurring in the host galaxy. 

In Sect. 2, we give a brief overview of the data and spectral 
models used, while in Sect. 3 we present the results. In Sect. 4, 
we discuss the obtained relative abundances, and compare these 
to the abundances quoted in the literature and derived using dif- 
ferent methods. Our conclusions are given in Sect. 5. 



2. Observations and spectral models 

2.1. Observations 

The data used in this study were obtained to determine the dis- 
tance from the ionising source of the absorbers in Mrk 509 by 
measuring the time delay between the change in flux and ioni- 
sation structure of the gas. This goal cons trained the observin g 
strategy, which is described in detail bv iKaastra et alj (1201 lal) . 
We present here a brief summary of the X-ray data used in this 
article, which consists of ten ~60 ks XMM-Newton observations 
taken four days apart on October-November 2009. The 180 ks 
LETGS (Chandra) spectrum was obtained 20 days later over two 
consecutive orbits between Dec. 10 and 13, 2009. 

Full details of the RGS da ta reduction, calibratio n, and stack- 
ing of the data are given in iKaastra et ^ (I2011bh and for the 
LETGS data by lEbrero et al.1 (1201 ll) . "V^note that the spectra 
were aligned using the very accurate aspect solution of the satel- 
lite before the stacking occurred, so that no line from the source 
is smeared out, and that special care was taken not to introduce 
weak emission or absorption features caused by a bad pixel in 
one of the spectra. 

2.2. Spectral models 

Here we summaris e the absorption model s given in 
Det mers et al.1 (1201 ih and lEbrero et afl (l201lh which we 
used to determine the relative abundances. The continuum 
was fit by a logarithmically spaced spline function, the ad- 
vantage being that this model does not make any a priori 
assumptions about the physical components of the continuum. 
Several narrow and broad emission lines were fit, as well as a 
few radiative recombination continua (RRC). Together these 
components allow for an accurate continuum fit, necessary to 
perform accurate line depth measurements. The local Galactic 
foreground absorption is fit with three collisionally ionised hot 
components with different temperatures. The hot model gives 
the transmission of a layer of collisionally ionised gas. Free 
parameters are the hydrogen column density, the abundances, 
and the temperature. 

The AGN absorption observed in the stacked RGS spectrum 
was fit with several different models. Model 1 includes two slab 
components, while model 2 uses three slab components. The 
slab model gives the transmission of a layer with arbitrary ionic 
composition. The main free parameters are the ionic column 
densities, but also the velocity outflow and velocity broadening 
are free parameters. In model 2 the ionic column density is deter- 
mined for all ions for the three velocity components. The ionic 
column densities obtained with model 1 were used as the input 
to the absorption measure distribution (AMD) modelling, which 
is the result we use in this paper. Model 3 fits six xabs compo- 
nents. The xabs model gives the transmission of a photoionised 
layer, where the ionic column densities are determined from the 
total hydrogen column density, the abundances, and the ionisa- 
tion parameter. For the LETGS, similar models for the absorber 
were used, although, for the xabs and AMD modelling only 
three components were detected significantly in the spectrum, 
owing to the lower signal-to-noise ratio (hereafter SNR). Owing 
to the complexity of the models used, one of which, model 2 for 
the RGS, has 22 components and nearly 100 free parameters, we 
decided not to fit the RGS and LETGS spectrum simultaneously. 
We instead fitted the abundances for the RGS and LETGS data 
separately, and compare in Sect j4.1l the values obtained for the 
separate analyses. 
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To obtain an acceptable fit to the RGS spectrum, we need 
at least three velocity components, and find consistent results 
for the LETGS spectrum. We note that for Mrk 509, seven 
differ ent velocity components are detected in the FUSE spec- 
trum (i Rriss etaUbOOO ). eight in the HST-STIS UV spectrum 
( Kraem er et al. [|2003h . and thirteen in the HST-COS spectrum 



SPEX softwar43 (iKaastra et alJI 19961) . and the errors quoted are 



(IKriss et al.l 201 ih . with velocities ranging between -426 and 



+219 km s' . The current X-ray spectra do not have the spec- 
tral resolution to separate these different velocity components. 
However, in the X-ray band we only detect absorption from ve- 
locities that are consistent with being at rest with the host galaxy 
or outflowing. The absorption component with the highest red- 
shift relative to the systemic velocity is only observed in the UV 
and is hkely to b e a high velocity cloud near the host galaxy 
dKriss et al.l 201 ll) . From studying the absorption troughs either 
at rest or outflowing in the UV, it appears that these are clustered 
around -300 and +16 k m s"' (iKriss et alj '2000'; 'Kraemer et al. 
120031: IKriss etal.ir201 lb . Fitting model 2. lOetmers et al.. (2011.) 
finds outflow velocities of -13+1 1 km s"' and -319+14 km s"^', 
consistent with the main absorption troughs detected in the UV. 
[Oetmers et al. (2011) also detects a higher velocity component 
at -770+109 km s"', but observed in only two ions: Mgxi and 
Fe XXI, that is not present in th e UV spectra. F i tting t he six differ- 
ent xabs models for the RGS. lDetmers et all (1201 lb indeed find 
that there is a slow (fit by components Bl and CI) and a faster 
(fit by components A2, C2, D2 and E2) compo nent, withdiffer- 
ent outflow velocities. In their model l. lEbrero et al.l ( l201ll) used 
two velocity components in the slab modelling similar to model 
1 for the RGS. Owing to a difference in absolute wavelength cali- 
bration between the RGS and the LETGS the two velocities are: 
-338 km s"' and +3 kms"', which are within the uncertainty 
consistent with the velocities fitted for the RGS and the two UV 
detected troughs. 

The stacked 600 ks high SNR RGS spectrum was analysed 
by iDetmers et al.l (l201 ll) . In this paper, we use the results from 
their model 3, which has six xabs compon ents and is detailed i n 
their Table 5, and AMD fit (see Table 6 in lDetmers e'tal]l201 ih . 
to determine the relative abundances. The latter is based on the 
results of model 1. For the LETGS, we use th e results of the 
xabs modelling obtained bv lEbrero et alJ (1201 ll) . as well as the 
AMD modelling. The only parameters we allowed to vary in our 
fits were the abundances of C, N, O (in the AMD fit), Ne, Mg, Si 
(only for the LETGS spectrum), S, Ca, and Fe, and the hydrogen 
column density. We thus kept the ionisation parameter, or the 
ionic column densities for the AMD fit, as well as the outflow 
velocity and th e velocity broa d ening fixed to the best-fi t value s 
determined by iDetmers etal] (120111) and lEbrero et al"] (1201 ll) . 
We did test whether leaving the ionisation and velocity parame- 
ters free significantly alters the derived relative abundance of S, 
but found that these were unchanged to one standard deviation, 
namely a change from 0.57 (see Table [T]) to 0.63, where the de- 
termined error in S abundance is 0.14 for the xabs modelling of 
the RGS. For the AMD modelling, we used the fitted ionic col- 
umn densities as input parameters, which hence have to be kept 
fixed. However, these ionic column densities were determined 
without assuming either an ionisation structure or abundances 
for the plasma. We briefly note that the continuum was signifi- 
cantly lower during the LETGS observation, 24.19 versus 37.78 
in 10 W m between 5 A and 40 A (iKaastra et alj|201 lal) . In 
addition to the shorter exposure time, this leads to a lower SNR 
in the LETGS spectrum. Therefore, the best-fit absorption and 
emission model is necessarily less detailed than that obtained 
for the RGS spectrum. All the spectral analysis is done with the 



Icr. 



3. Derived abundances 

Hydrogen only produces continuum absorption, which is de- 
generate with our continuum parameters. Thus determining the 
absolute abundances, i.e. the abundances relative to hydrogen, 
which would likely be in the form of H n, is difficult. Hence, we 
used abundance ratios where we chose an element with a large 
number of absorption lines over a large ionisation range and de- 
termined the abundances relative to this element. For this ele- 
ment, we used oxygen because it is the most abundant metal and 
we observe transitions from O iv to O vm with multiple detected 
lines per ion, with the exception of O iv. 

3. 1 . RGS determined abundances 

3.1 .1 . Abundances determined assuming all components 
have the same abundances 

The stacked RGS spectrum has a much higher SNR than LETGS 
and thus provides us with more accurate relative abundances. We 
use model 3, which fits the absorption with six xabs compo- 
nents that each have a different outflow velocity and the AMD 
modelling to determine the relative abundances. For the parame- 
ters used in these models we refer to IDetmers et al. (201 I p, who 
assumed the proto-solar abundances given by Lodders & Palmd 
(I2OO9I) . Here we determine the relative abundances using three 
methods based on model 3, the AMD model, and different as- 
sumptions. 

We first use model 3 and assume that all the differently 
ionised and different kinematic components have the same abun- 
dances. We allow the abundances of C, N, Ne, Mg, S, Ca, and 
Fe to vary and choose O as the reference element. Since no Si 
absorption line is visually detected in the RGS spectrum, for 
the RGS abundance determination we did not fit the Si abun- 
dance. We thus couple the abundances for the six xabs compo- 
nents. The resulting relative abundances are given in column 2 
of Table [T] From Table [T] we conclude that the abundance ratios 
are consistent with the proto-solar ratios with the exception of 
the S/O. 

For iron, which has a best-fit value 2.5cr below the proto- 
solar ratio, some of the stronger lines are notably weaker in the 
best-fit model, but the main difference is in the absorption due 
to the unresolved transition array (U TA between 16 and 17.5 A, 
ISako et aTl 120011: iBehar et al]|2001b . because the many absorp- 
tion lines result in pseudo-continuum absorption. The spacing of 
2.68 A and 3. 18 A between the points of the spline used to fit the 
continuum in the range 14-18 A is wider than the UTA, which 
has a width of about 1.5 A. 



3.1 .2. Abundances of each velocity component 

Thus far we have assumed that the faster and slower velocity 
components and the different ionisation components have the 
same abundance. The latter assumption seems justified, if the 
different ionisation components are part of the same outflow and 
thus form one physical entity with a common origin. The former 
assumption is more questionable if the different velocity com- 
ponents are from different outflow structures formed at different 
distances and possibly with a different origin. The assumption 
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Table 1. Best fit abundances relative to oxygen compared to the proto-solar abundance ratio from the RGS spectrum (using model 
3 and the AMD models, see Sect l3.1b and the LETGS spectrum using the xabs and discrete AMD model (see Sect l3.2b . We assume 
that the abundances for the faster and slower outflow component are the same. 



ion 


RGS: 
model 3 


discrete AMD 


continuous AMD 


LETGS: 
xabs 


discrete AMD" 


recommended* 


ion' 


ion'' 


C/O 


1.19±0.08 


1.00+0.07 


1.01+0.03 


1.2+0.3 


1.00 


1.19+0.08 


C/Fe 


1.40+0.1 


C 


8.46 


N/O 


0.98±0.08 


0.80+0.05 


0.80+0.03 


0.5+0.2 


0.80 


0.98+0.08 


N/Fe 


1.15+0.06 


N 


7.90 







1.00+0.05 


1.00+0.05" 




1.11 




O/Fe 


1.17+0.06 





8.76 


Ne/O 


1.11+0.10 


1.17+0.15 


1.16+0.07 


0.9+0.3 


1.17 


1.11+0.10 


Ne/Fe 


1.31+0.17 


Ne 


7.95 


Mg/O 


0.68+0.16 


0.92+0.29 


0.93+0.15 


0.7+0.6 


2.9+5.3 


0.68+0.16 


Mg/Fe 


0.80+0.16 


Mg 


7.62 


Si/O 








1.3+0.6 


1.2+1.5 


1.3+0.6 


Si/Fe 


1.53+0.60 


Si 


7.61 


S/O 


0.57+0.14 


0.68+0.24 


0.69+0.09 


0.4+0.4 


1 + 1.1 


0.57+0.14 


S/Fe 


0.67+0.15 


S 


7.26 


Ca/0 


0.89+0.25 


3.34+0.75 


3.81+0.86 


2.7+1.6 


5.1+4.7 


0.89+0.25 


Ca/Fe 


1.05+0.26 


Ca 


6.41 


Fe/O 


0.85+0.06 


0.81+0.05 


0.79+0.08 


1.1+0.2 


0.81 


0.85+0.06 


Fe 




Fe 


7.54 



" Values without error bars denote that they were kept fixed to the best fit discrete AMD model of the RGS spectrum. 
* The most accurately determined relative abundances that we use in the discussion of this paper. 
'■ The calculated relative abundances relativ e to Fe, for com parison with other measurements. 



'' Logarithmic proto-solar abundances of Lodders & Palme (2009), where H has a value of 12. 
" Fixed to the value obtained for the discrete AMD model. 



that all o utflow component s have the same abundance was also 
made by lArav et al] (l2007h in their study of the abundances of 
Mrk 279 using FUSE and HST-STIS UV spectra. The lower ve- 
locity outflow component might be contaminated by absorption 
from the host galaxy. However, as we observe the host galaxy 
face on and detect no absorption from neutral gas in the X-ray 
spectra at the redshift of Mrk 509, this should have a limited 
effect on the derived abundances. However, we note that in the 
UV substantial C m and Siiv is detected from the host galaxy 
jKriss et al. and there could thus be contamination, cer- 

tainly for the more lowly ionised gas. Another possib le contami- 
nant is the inflowing gas seen in UV absorption lines (iKriss et alJ 
12011 ). We note that in the UV detected Ovi line the strongest 
component has a velocity of -1-219 km s"' but this velocity com- 
ponent is not detected in the RGS spectrum. Therefore, the con- 
tamination by absorption from inflowing gas is also likely to be 
limited. 

Considering the small error bars in the determined abun- 
dances, we decided to attempt, as a second method the determi- 
nation of the abundances separately for the fast and slow outflow 
components. We still tied the abundances for the different ioni- 
sation components for each velocity component. We again use O 
as a reference element. We thus fit two abundances per element 
simultaneously. As expected, the error bars are larger owing to 
the correlations between both velocity components. 

With the exception of the Fe/O ratio for the slow velocity 
component, all abundance ratios are consistent with proto-solar 
ratios and the values determined assuming the abundance is the 
same for both velocity components. The only difference is for 
iron, but this might be due to the difficulty in fitting the many 
weak and lowly ionised lines, as well as some uncertainty in the 
centroids of some of the blends. 

3.1.3. Abundances from the absorption-measure distribution 
models 

Since the effective spectral resolution is insufficient to separate 
the slow and faster velocity components into independent com- 
ponents, and that likely part of the fastest outflow velocity de- 
tected in just two ions is improperly fit with the six xabs com- 
ponents, we decided to use the ionic co lumn densities deter- 
mined from model 1 dDetmers et al.ll20Tll) to derive the absorp- 



tion measure distribution and the abundances simultaneously. 
The procedure used to determine the AMD taking into account 
the full ionisation range at which an ion is formed is described 
by Detmers et al. (2011). 

We used in this method, in contrast to the previous methods, 
the ionic column densities measured from the spectrum using 
two slab components. Thus, a priori we did not assume any 
ionisation structure, which in model 3 is assumed to be dis- 
crete and determined by the xabs components. To obtain the 
relative abundances from the measured ionic column densities, 
we did need to assume a certain SED and ionisation balance, 
which were taken to be the same as for the xabs modelling. 
An advantage of this method is that we can test whether the 
ionisation str ucture is discrete or cont inuous, as suggested for 
NGC 5548 bv lSteenbru^ge et al.l(l2005h . We thus determined the 
abundances for both discrete and continuous ionisation models. 
The details of these fits are described in Detmers et al. (201 1|. 
Since no hydrogen lines are present, we were unable to deter- 
mine absolute abundances. We determined the abundances of 
individual elements relative to all other elements with measured 
column densities, because the fit first derives the total column 
density assuming proto-solar abundance ratios and the fit is then 
refined by adjusting the abundances. As before, we therefore re- 
normalise our abundances to oxygen. The error bar given for 
oxygen is the nominal uncertainty in the oxygen abundance rel- 
ative to all other elements. 

The abundances for a discrete and a continuous AMD model 
are given in columns 3 and 4, respectively, of Table[T] Consistent 
abundances were obtained by assuming either a discrete or a 
continuous ionisation structure. This provided us with confi- 
dence that we can determine the ionisation structure of the gas 
accurately e nough to prevent any i nfluence on our abundance de- 
termination. iDetmers et al.l (1201 lb does convincingly show that 
the ionisation structure is discrete, at least for highly ionised gas, 
hence we only discuss this model further here. 

3.2. LETGS determined abundances 

Here we repeated two of the above mentioned methods to mea- 
sure the abundance, but this time for the Chandra LETGS data. 
Owing to the lower SNR of the data, we used C-statistics to fit 
the data, instead of statistics as was the case for the RGS 
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data. We determined the abundances usin g the spectral mode l 
with three xabs components described by Ebrero et"an (1201 ih 
assuming the continuum is well described by a spline. The dif- 
ference in absorption parameters are negligible if a power-law 
and modified black body were assumed instead of a spline. The 
outflow velocity and the velocity width were fitted separately for 
each xabs component, as done for the RGS analysis. Owing to 
the lower SNR, only three instead of six xabs components are 
significantly detected (Ebrero et al. 2011 ). This is in part because 
the slower and faster outflow components are only partially re- 
solved and because the lowest ionisation component is not sig- 
nificantly detected. The lowest ionisation component also has 
the lowest hydrogen column density, and is thus more difficult 
to detect significantly in lower SNR data. We do know that low 
ionisation gas was present during the LETG S observation, be - 
cause it was detected in the HST-COS spectra (iKriss et al .1201 Ih . 
The best-fit abundances that we derived using three xabs com- 
ponents, and tying the abundances for the diff'erent xabs compo- 
nents, for the LETGS spectrum are listed in column 5 of Table[T] 

We also determined the abundances using the AMD mod- 
elling, assuming a discrete ionisation structure, using the ionic 
column densities determined for mo del 1, which includes two 
slab components (lEbrero et al l20Tlh . This model is the same as 
model 1 for the RGS and has two different velocity components, 
although each ion is only fit using one velocity component. The 
AMD fit only detects three discrete ionisation components, ow- 
ing once again to the lower SNR of the data set, and consistent 
with the xabs modelling. Considering the lower SNR, we as- 
sumed that the ionisation structure of the gas is discrete, as this 
gave a more accurate description of the RGS data than the con- 
tinuous model. However, the fit is unstable if we allow the ion- 
isation parameters, hydrogen column densities, and abundances 
to be free parameters, owing to the strong correlations between 
these parameters caused by the limited SNR. The error bars in 
the ionic column densities are too large, thus we decided to fix 
the abundances for those elements that have a well-determined 
abundance from the AMD analysis of the RGS data (C, N, O, 
Ne, and Fe). However, owing to the larger wavelength range, we 
can determine the Si abundance, and possibly more tightly con- 
strain the Mg, S, and Ca abundances. 

The uncertainties in both methods are generally quite large 
for the ions for which we allowed the abundance to vary, be- 
cause these elements show only a few weaker lines and the abun- 
dances were also poorly determined with the RGS spectrum. The 
abundances of these elements (Mg, Si, S, and Ca) are consistent 
with proto-solar abundances. The improvement in the fit with the 
best-fit abundances and the proto-solar abundances using xabs 
modelling is small, indicating that they are not very secure de- 
terminations, as can also be seen from the error bars and the dif- 
ference in the values derived by the xabs and AMD modelling. 
As a result of the larger error bars, the relative abundances are 
consistent with those determined for the RGS spectrum. 

Since the higher SNR and the spectral resolution of the RGS 
spectrum are too poor to accurately determine the abundances 
for each velocity component, we did not try to fit the abundances 
for each velocity component in the LETGS spectrum. 

4. Discussion 

We used the narrow absorption lines to measure the relative 
abundances, but instead of using the UV-detected narrow absorp- 
tion lines, as Arav et al. (2007) did, we used the X-ray detected 
narrow absorption fines. An advantage of using absorption lines 
is that the lines, with the exception of very weak and in the case 



of Mrk 509 undetected metastable lines, are independent of the 
density and temperature of the gas. 

The use of X-ray, instead of UV, detected absorption lines 
has several advantages. One advantage is that there are absorp- 
tion lines from enough ions to accurately determine the ioni- 
sation structure of the absorbing gas. The many iron and oxy- 
gen ions detected in X-ray spectra enable the ionisation struc- 
ture to be accurately determined using only the oxygen or iron 
ions. This can be done without abundance differences compli- 
cating the analysis, something that is impossible from the lines 
detected in the UV spectrum. The analysis of UV spectra is also 
complicated by the saturation of the absorption lines and the fact 
that one has to correct for non-black satura tion and a velocity- 
dependent covering factor dArav et al ] l2003h . In the X-ray spec- 
trum of Mrk 509, and for most of the other Seyfert 1 galax ies, the 
covering factor is unity (Detmers et al. 2011; Steenbrugg e et al.l 
i2009h . We did test whether the covering factor is indeed unity by 
forcing the covering factor of the six xabs components to be 0.7 
and refitting the RGS spectrum, allowing the hydrogen column 
densities and ionisation parameters of the different xabs compo- 
nents to vary as free parameters. The resultant increased by 
55 from the best fit with unity covering factor If we then let also 
the covering factors free, for all six components this resulted in a 
best fit with unity covering factors. Furthermore, in X-ray spec- 
tra there are enough line series and weaker lines to correctly de- 
termine the column density of the ions, even if the deepest lines 
are saturated. We should note that in Mrk 509 the lack of detec- 
tion o f the weaker absorption li nes observed in for instance NGC 
5548 (ISteenbrugge et al.]l2005h indicates that generally, i.e. with 
the exception of the strongest lines, the absorption lines detected 
are not saturated. A final advantage is the large range in ionisa- 
tion parameter sampled in the X-ray band, allowing us to study 
many ions of an element when determining its abundance. 

However, there are also disadvantages. One is that we cannot 
separate the velocity components detected in UV spectra, thus as 
in abundance determinations based on UV spectra, we have to 
assume that all velocity components have the same abundances. 
Another disadvantage is that in the X-ray band no transitions of 
H are observable, thus the determined abundances are all rela- 
tive. Only FUSE enabled the full Lyman series to be observed 
for this redshift range, permitting absolute abundances of AGN 
to be determined. 

The measured abundances relative to oxygen are consistent 
with proto-solar abundance ratios with the exception of sulphur 
For S, the different methods and the different instruments lead to 
a result implying that S is more likely to be under-abundant than 
have a proto-solar abundance ratio. For this element, we did test 
whether allowing the ionisation and velocity parameters to vary 
in the fit that determines the abundances has a significant effect 
on the determined S abundance ratio. However, the additional 
free parameters increase the abundance by less than 1 cr for the 
RGS spectrum fit with the six xabs components. The abundance 
is determined from three lines longward of 3 1 A, and it is thus 
a fairly robust result. Sulphur is formed mostly in supernovae 
(hereafter SNe) type la, as are Ca and Fe, although SNe II are a 
far more significant contributor to the S abundance than to either 
the Ca and Fe abundance. If we study the abundance ratios rel- 
ative to iron rather than oxygen, then S has an abundance ratio 
consistent with proto-solar, but carbon is slightly overabundant. 
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4.1. Comparison between the abundances derived 



4.2.1. Abundances from narrow absorption lines 



All the abundance ratios, apart from that of Ca, derived for 
the RGS spectrum are consistent between the different methods 
used. The differences between the relative abundances for the 
continuous versus discrete AMD and the abundances derived us- 
ing model 3, provide a good indication of the systematic errors 
caused by the different ionisation structure modelling. Since the 
abundance ratios inferred from the AMD modelling are consis- 
tent with those derived using model 3, owing to the error bars 
being on average larger, for the remainder of the paper we use 
the abundances determined from model 3, assuming that all the 
velocity components have the same abundances. 

The abundance measurements are consistent to within 2.5cr, 
and there is no general trend observed for the results of the RGS 
and LETGS data sets. The ions fit for both the RGS and LETGS 
data using the AMD model give abundances that are consistent 
to within Icr. This is unsurprising, as we only fit those ions that 
had larger error bars in the RGS abundance determination, in the 
LETGS spectral fitting. 

The Ca/O abundance measurements are slightly inconsistent 
to within 3cr between the different methods used to determine the 
abundance from the RGS spectrum. The xabs measured abun- 
dance is consistent with the proto-solar Ca/O ratio, while the 
measured AMD abundance is inconsistent with the proto-solar 
ratio. In both the RGS and LETGS spectra, only a Caxiv ab- 
sorption line is visually detected. This is also the only ion that is 
detected in the slab modelling (models 1 and 2) of the absorber. 
The detected Ca xiv line is blended with absorption from the lo- 
cal Nvii Lya line. Owing to this blending and since no Caxiii 
or Caxv lines are detected, and the xabs modelling is more re- 
liable in determining the column densities of weak and blended 
features, we use the RGS data xabs abundance in the rest of this 
paper. 



4.2. Comparison between our abundances and nuclear 
abundances in tine literature 



The chemical abundances of quasars provides insight into the 
evolution of their host galaxy and in particular the enrich- 
ment processes such as star formation. Therefore, the abun- 
dances of high redshift quasars especially have been well- 
studied. The two principal methods for measuring quasar abun- 
dances are the study of the broad emission Hnes (Shields 19761 
[Hamann & Ferland 1999) a nd the study of the narrow absorption 
hnes (iGabel et al.i,2006i: lArav et al.ll2007h . The broad emission 
Hnes are formed in the vicinity of the supermassive black hole at 
distances from the black hol e generally between 1 -70 light-days 
(8.4x10^^ - 0.059 pc) ( Peterson & Wandel 2000). 

The outflow also originates close to the super-massive black 
hole, as it is released either by the accretion disk or the torus. 
Hence, both methods determine the abundances in the vicinity 
of the central black hole. Between redshifts two and four, the 
study of the broad emission lines infer N/ C ratios of up to 2, 
hence metallicities of up to 14 times solar dHamann & Ferlaiidl 
[T99I . Even for local quasars (z<0. 1), the av erage metallicity de- 
rived from the N/C ratio is about twice solar (Hamann & Ferland 
Il992h . Using the second method, Arav et al. (2007) measured 
absolute abundances for Mrk 279, which they found to be con- 
sistent with either solar or twice solar values. 



Narrow absorption lines we re used to d e termin e the abundances 
of the observed outflow by Ara v et aTl ( l2007h in the UV spec- 
trum of Mrk 279. These authors measured the absolute abun- 
dances of C, N, and O to be C 2.2+0.7, N 3.5+1.1, and O 
1.6+0.8, corre sponding to 1.9+0.7, 3.6+1.1, and 1.3+0.8 on the 
iLodders & Pa lme (2009) scale. If we compare our relative abun- 
dances to the absolute abundances measured for Mrk 279, we 
find that the C/O and N/O ratios are somewhat larger for Mrk 279 
than for Mrk 509. However, their error bars ar e also larger, thus 
the abundance ratios measured for Mrk 279 dArav et alj|2"00l 
are consistent within the error bars with both proto-solar val- 
ues and the abundance ratios measured for Mrk 509. Using the 
VLT UVES spect rum of the AGN w ind of the z=2.238 quasar 
QSO J2233-606. iGabel et alJ (l2006h measured super-solar ab- 
solute abundances of C 4.0-5.0, N 12.6-15.9, and O 3.2-5.0. 
The abundance ratios thus range between 1 and 1 .6 for C/O and 
2.5 and 5 for N/O, indicating that there is a clear overabundance 
of N relative to O, not observed for Mrk 509. 

It would thus appear from the abundance measurements in 
Mrk 279 and Mrk 509 derived using high-resolution spectra, that 
in the local universe, the abundance ratios of the gas near the 
accretion disk are close to proto-solar/solar values. 

There have been several previous attempts to determine rel- 
ative abundances using X-ray spectra, but generally the SNR 
of the data were too low for these measurements to be con- 
strained. For earlier results, inaccuracies in the atomic data used 
also led to biased results, especially for Fe when the most re- 
cent di-electronic recombination rates w ere not used, as dis- 
cussed by Netzer (2004). For instance, Sc hmidt et alj (l2006h and 
INahar et alj ((2009) performed detail ed atomic cal culations of 
the di-electronic recombination rates. .Schmidt et an(2004 ) mea- 
sured the wavelengths of the O v Unes. after ISteenbrug ge et al.l 
(120031) had noted that the wavelength used for the 22 A line was 
likely to be wrong. In their appendix 1, Detmers et al. (201]]) 
list the atomic data that has been updated prior to the analy- 
sis of the RGS and LETGS spe ctra of Mrk 5 09. For the 40 ks 
RGS spectrum of NGC 3783, iBlustin et al] (|B002) measured 
the abundances to be consistent with solar, with th e exception 
of Fe. Similar abundances were measured by Steenbrugge e t al.l 
( I2003h for NGC 55 48 using a 137 ks RGS spectrumras well 
as by IBlustin et all (l2007i) for a 160 ks RGS spectrum of 
NGC 7469. iHolczer et alJ (l2007h determined the abundances 
for IRAS 13349+2438 and NGC 3783 using the updated di- 
electronic recombination rates and determined abundances rel- 
ative to Fe consistent with proto-solar ratios for O, Ne, Mg, Si, 
and S (NGC 3783 only). In Table |2l we compare our relative 
abundances to those obtained by Holczeret al. (2007), convert- 
ing the abundances to the iLodders & Paknd ([2009) abundance 
scale. Although some of these earlier abundance measurements 
suffer from large uncertainties, the determined abundances of 
AGN in the nearby universe using narrow X-ray absorption lines 
are broadly consistent with solar abundances. 

4.2.2. Abundances from broad emission lines 

We can also compare the abundances we derived, and those de- 
rived for the outflow in other nearby Seyfert 1 galaxies, with 
those derived for the broad emission-line region. Although the 
relationship between the absorbing gas and the BLR is unknown, 
since both are located near the accretion disk, their abundances 
should be similar, and we proceed under that assumption. In the 
approach that we adopt, some additional assumptions are nec- 
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Table 2. Comparison of the relative abundances of Mrk 509 
(see Table [T} with those derived for IRAS 13349+2438 and 
NGC 3783 (Holczer et al. 2007). 



Mrk 509 IRAS 13349+2438 NGC 3783 



O/Fe 


1.17+0.06 


1.3+0.4 


1.0+0.2 


Ne/Fe 


1.31+0.17 


1 3+"' 

i-'_0.4 


1 4+0.3 


Mg/Fe 


0.80+0.16 


1.4+0.5 


1 5+°" 
li^:^ 


Si/Fe 


1.53+0.60 


2.2+0.75 


S/Fe 


0.67+0.15 




1 0+^-^ 



essary to determine the abundance, as in the optical and UV 
spectra only a handful of broad emission lines are observed. 
As a result, the ionisation structure of the gas is poorly con- 
strained, and one thus has to rely on a model for the ionisation 
structure and to use lines with the same ionisation parameter, 
otherwise there is a degeneracy between the ionisation struc- 
ture and abundance measured. In these studies, it is generally 
assumed that the N/O ratio is a proxy for O/H, which is itself 
a proxy for the met allicity Z, or that N/H is p roportional to the 
metallicity squared jHamann & FerlandllI999h . The relation be- 
tween N/O and meta llicity was derived and is verified for Hn 
region s (ISearlelll971h . It was applied to quasar nuclei bv^ Shieldsl 
( Il976l) . who also expanded the relationship by indi cating that 
N/C should follow the same trend as N/O. Shields( (1 19761) and 
other early observers found a significant nitrogen overabundance 
relative to oxygen, using the N iv] and N m] lines detected for a 
handful of nearby active galaxies; however, for the average spec- 
trum of quasars they did m easure N/O ratios of one-fourth to half 
solar. iHamann & Ferlandl Cl993) showed that the N v/Civ and 
the N v/He n ratios are less dependent on the ionisation parame- 
ter and are thus to be preferred. The N v, C iv, and He ii lines are 
therefore currently the most widely used lines in determining the 
abundances from the broad emission lines. 

Using the N v/C iv and N v/He n ratios, high redshift quasars 
are found to have a range of abundances between a few 
to ten times solar (.Ham ann & Ferland 1999). These authors 
also derived a sfi ght in crease in abundance with redshift. 
IHamann & Ferlandl dl992l I1993L [T999.) conclude therefore that 
quasar activity only becomes detectable after vigorous star for- 
mation, although the enrichment could be so rapid that star for- 
mation and QSO formation are coeval. This is in contrast to the 
earlier results, based on diff'erent emission lines and the study of 
lower redshift quasars, that generally found quasar abundances 
are sub-solar Since both methods use gas that is formed at a very 
similar distance from the inner accretion disk that is assumed to 
have the same origin as the disk, the difference in abundance is 
due to either the assumed ionisation structure for the gas, be- 
cause of galaxy evolution, or because one, or more, of the lines 
used in the abundance calculation is contaminated. Below we 
describe the observational evidence against such high metallic- 
ities at high redshifts, and argue that the N v fine used is most 
likely contaminated. 

The high metallicities measured at high redshift imply 
that the initial mass function (IMF) was relatively flat dur- 
ing an intense period of star formation, which thus pro- 
duced mor e high-mass stars than ob s erved in the solar neigh- 
bourhood (IHamann & Ferlandl Il999l Il993h . This is inconsis- 
tent with the IMF derived for the eight most luminous and 
massive cluster galaxies in t he nearby universe studied by 
Ivan Dokkum & Conrovl (l2010l) . These authors derive a steeper 
than Salpeter IMF, namely a power-law slope of — 3. This is 



much steeper than the power-law slope derived for the Milky 
Way disk, and implies that 80% of stars in these galaxies are 
low-mass stars, van Dokkum & Conrov (2010) do point out that 
this IMF is applicable also to the central region, where the ini- 
tial starburst should have occurred and that star formation caused 
by later mergers mostly affect the stellar mass function at larger 
radii. It thus seems that the high metallicities measured from 
the broad line region are inconsistent with the observed IMF in 
bright massive nearby cluster galaxies. 

The possibility that t he N v line is contaminated was first 
suggested by Turnshekj (Il988l) . who found a possible corre- 
lation between the N v broad emission line strength and the 
broad absorption li nes measured in the AGN under study. 
iKrolik & Voill (Il998l) calculated the possible contribution of res- 
onant scattering t o the N v line to be 0.4-0.5 times the Lya 
flux. Finally, Wan g et all (I20T0I) summarised the observational 
evidence which implies that indeed the Nv line is generally 
strongly contaminated by the resonance scattering of Lya pho- 
tons. Strong support is p rovided by the 10% polarisati on fraction 
of the N V emission line ("Lamv & Hutsemekersll2004l) . while the 
finding that other broad emission lines are unpolarised indicates 
that at least part of the emission is scattered light. Furthermore, 
IWevmann et al.l ( 1199 lb found a difference between the strength 
of the Nv line in BAL and non-BALQSOs, with excess Nv 
emission extending to 10,000 km s"', much broader than the 
broad emission lines. The contamination with resonant-scattered 
Lya photons diminishes the reliability of the metallicities de- 
rived from the broad emission line region using the N v line and 
might lead to spurious over-abundances. However, when using 
other line ratios the uncertainties caused by the difference in the 
line ratios with temperature, ionisation structure, density, and as- 
sumed spectral energy distribution ensure that the abundance de- 
termination is uncertain. Another complication of the abundance 
determination from broad emission lines is that the BLR is strati- 
fied, i.e. the ionisation parameter and hence line strengths change 
with distance from the black hole. Nevertheless, the abundances 
derived for the BLR are overwhelmingly using one zone models. 
Therefore, we conclude that the abundances determined for the 
BLR are unreliable, and especially when the N v line is used. 

4.3. Abundances in the Galactic centre 

Our own Galaxy is the only case where we can distinguish the 
abundances of the different components in the vicinity of the su- 
permassive black hole. As a result, the abundances have been 
determined for both molecular clouds and H n regions and from 
both cool and hot stars. For the molecular clouds in the vicinity 
of the Galactic centre (GC), the derived iron abundance depends 
on the assumed process produ cing the observed X-ray emis- 
sion, and ranges from 1.3+0.1 (iTerrier et alJl2010l) . to 1.6+0.2 
(Revnivtsev et al.''2004'), and to 3.4 times the pro t o- solar value 
( Yusef-Zadeh et al. 2007). Studying Hn regions. fOiveon et al.l 
(2002) determined a metallicity gradient for the Galaxy, which 
gives an abundance of the GC for Ne of 2.3+1.3 and for S of 
0.68+0.27 times the proto-solar value. The Ne/S ratio of the H 11 
regions near the GC is larger than that measured for Mrk 509, but 
both measurements are consistent. From the stud y of 1 1 gi ants 
and supergiants within 30 pc of the GC, ICunha et al.i (l2007t) de- 
rived an average proto-solar iron abundance of 1 . 12+0. 17 and an 
average O abundance of 1.9+1.1. Solar zero-age main-sequence 
abundances were derived from the study of five late-type Wolf- 
Rayet stars by Naiarro et al. (2004). Hence, most of the measure- 
ments of the absolute abundances in the immediate surroundings 
of Sgr A* are consistent with proto-solar and twice proto-solar 
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abundances. The O/Fe abundance ratio derived from cool giant 
and supergiant stars is consistent with the measured abundance 
ratio of Mrk 509. 



4.4. Metallicity of star-forming galaxies 

Here we compare the abundances measured for Mrk 509, with 
measurements of the metalHcity in the interstellar medium (ISM) 
in nearby star-forming galaxies. We note that Mrk 279 is the 
only Seyfert 1 galaxy with reliable absolute abundances for the 
outflow. The gas near the accretion disk is due to inflow, hence 
the measured metallicity should be related to the metallicity 
of the ISM. In the comparison, we implicitly assume that the 
abundances of the gas falling toward the black hole do not dif- 
fer significantly from the average abundances/metallicity of the 
host galaxy. Ideally, we would compare the measured nuclear 
abundances with the ISM abundances of the host galaxy of the 
Seyfert galaxies under study, but the metallicities of these galax- 
ies are unavailable. Nor is there a known mass-metallicity re- 
lationship for either Seyfert or active galaxies. There fore we 
use th e mass-metallicity relationship derived by Tre monti et al.l 
(l200l for star-forming galaxies, in addition assuming that the 
star-forming galaxies in the local universe have had a simi- 
lar enrichment history to the host g alaxies of Mrk 279 a nd 
Mrk 509. The metallicity is defined in lTremonti et al.l ([2004') as 
12-i-log(0/H), where a value of 8.69 is the solar oxygen abun- 
dance d etermined by ^ llende Prieto et alj (l200l . To compare 
with the iLodders & Palmg (i200p ) abundances u sed, in the fol- 
lowing we divide the values obtained using the iTremonti et alJ 
( |2004|) relationship by 1.235. We do not have absolute abun- 
dances for Mrk 509, making a direct comparison i mpossible. 
Howe ver, for Mrk 279 the absolute O abundance (lArav et alJ 
l2007 h is available and a comparison is possible. 

Mrk 509 is a very compact galaxy and therefore we assume 
that the host galaxy consists of only a bulge. We can calculate 
th e mass of the gala x y to be ~5x 10'° M© using the relation given 
by iMagorrian et al.l (Il998h an d using the measured m ass of the 
black hole Mbh = 3x10** Mq dMehdip our et al."2011). The un- 
certainty in the mass of the black hole obt ained from reverbera- 
tion mapping is about 30% dPeterson et al.ll2004l) . and is higher 
for the mass of the whole galaxy owing to the spread in the 
Magorrian relation as well as our assumption that t here is only 
a bulg e in this clearly disturbed galaxy (see Fig. 8 of lKriss et alj 
1201 Ih . However, the metallicity determination is only slightly 
affected by assuming masses of only 2x10'" or 8x10'" M©. 

For a galaxy mass of 5x10'" M©, we obtain an O/H ratio that 
is 1.68 times the proto-solar value. For a mass of 2x10'" Mq the 
O/H ratio is 1.65, and for a mass of 8x10'" Mq it is 1.70 times 
the proto-solar value. The uncertainty in the assumed galaxy 
mass is rather large, but cannot explain the difference between 
the nuclear O abundance and the higher global abundance de- 
temined for this galaxy. There is a spread in metallicity along 
the assumed relation of 0.1 dex, which is also unable to ex- 
plain the supe r -proto -solar abundance. The fiber used in the 
iTremonti et alJ (l2004l) study does not cover the whole galaxy, 
hence the integrated metallicity for the whole galaxy should be 
somewhat lower. The method of using nebular lines to deter- 
mine m^talli^itYjmgh^ the metallicity by a factor of 
two dKennicutt et al.ll2b03h . If this were indeed the case then the 
metallicity of the host galaxy ISM would be about proto-solar. 
Determining absolute abundances and the relative abundances 
of different elements is difficult using nebular lines, for the same 
reasons as the measurements of abundances from broad emis- 
sion lines is, which include the uncertainties in ionisation, tem- 



perature, and density of the gas, and that an unknown fraction of 
the heavy elements is locked up in dust. Since we only measure 
abundance ratios relative to O, the derived metallicity cannot be 
directly compared with our measured abundances for Mrk 509. 
However, we do find abundance ratios of the different elements 
are close to proto-solar values. 

lArav et alJ (l2007h did use the Allen de Prieto et al.1 (1200 ih 
abundances, thus for Mrk 279 the metallicity could be directly 
compared to the absolute abundance. Th e mass of the host 
galaxy bulge of Mrk 279 is 2.29x10'" Mq ( Wandel 200 2i), and 
the g alaxy is classified as a SO galaxy (de Vaucouleurse t al.l 



1199 Ih . Assuming that the mass of the disk is 30% of the bulge 
mass, the galaxy mass is somewhat lower than the bulge mass 
calculated for Mrk 509, ~3xlO'" Mq. Hence, the O/H ratio is 
1 .66, not corre cting for a p ossible factor of two overprediction of 
the O/H ratio. lArav et all ^bO[7j measured a O/H abundance of 
1.3+0.8, which is consistent with either proto-solar or twice the 
proto-solar metallicity. The large uncertainty in both the mea- 
sured abundance and the galaxy mass determination, as well as 
the possible factor of two overprediction in metallicity, imply 
that the results are inconclusive, and clearly more reliable data 
are needed to improve this comparison. 



4.5. Abundances from tfie fiot ISM 

An alternative way to study ISM abundances is to use X-ray 
emission lines from the hot ISM observed in some local ellip- 
tical galaxies. Owing to the larger number of ions and lines that 
can be observed, the problems in determining ionisation, temper- 
ature, and gas density are less severe than in the optical for the 
broad emission lines. However, these abundance determinations 
depend on the assumed thermal structure, as it is well known that 
a single temperature model p redicts t oo low abundances if the 
gas has two temperatures (iTrinchieri et al.lll994t iBuotdbOOOl) . 
Furthermore, in a few cases resonance scattering needs to be 
taken into account for the strongest iron lines. At the outskirts of 
the galactic X-ray halo, there might be contamination by either 
cluster or group emission, although, this does not affecUhe abun- 
dances derived for the centres of these haloes. |ji et al.l (l2009h 
studied ten local elliptical galaxies that have a bright X-ray halo 
to determine the abundances of the hot interstellar medium. The 
median relative abundances, for the elements that we also mea- 



sure in Mrk 509, in the eight galaxies for wh ich iJi et al.l (120091) 
obtained an acceptable fit are (converting to ILodde rs & Paling 
(2009) abundances) O/Fe = 0.66, Mg/Fe = 0.94, and Si/Fe = 
1.6+0.6. These should be compared with our measurements of 
1.17+0.06, 0.80+0.15, and 0.67+0.16, respectively, where the 
Si abundance is taken from the xabs modelling of the LETGS 
spectrum. Although for the individual galaxies, Iji et al ] (I2009h 
determined a wide range of abundances there is a rather large 
difference in the O/Fe abundance compared to Mrk 509, for 
which the 75% quartile is 0.77. The Mg/Fe and Si/Fe abundances 
are consistent with our results for Mrk 509. The oxygen abun- 
dance is generally observed to be low compared to the Fe and 
Mg abundances in the hot interstellar medium (Iji et alJ [20091 : 
Humphrev & Buote 2006). The reason for this low O abun- 
dance is a bit of a mystery, but might indicate that massive stars 
lost more of their He-rich atmosphere before the supernova ex- 
plosion (hereafter SN) occurred, and hence that there was less 
He burning in the outer envelope than current theories predict 
(Fulbright et al. 2005). For two galaxies of the ten studied by 
IJi et al.r ( f2009t) . the O/Fe ratio determined with the RGS in the 
inner 1' is close to the proto-solar ratio, hence our results are not 
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necessarily inconsistent with the abundances measured for the 
hot interstellar medium of elliptical galaxies. 

4.6. Comparison with cluster abundances 

Finally we also compare the abundances to those for the intra- 
cluster medium (ICM) in clusters of galaxies. The large gravita- 
tional potential of the cluster allows it to retain the gas ejected 
from the galaxies forming the cluster. It has generally been as- 
sumed that from the intra-cluster abundances the relative ra- 
tio of SN type la versus core-collapse SN can be constrained; 
however, stellar winds do play a role (Simion escu et alj|2009b . 
There has been some debate about the mixing of the differ- 
ently enriched ejecta, but most metals generally appear to have 
a peak abundance in the core, even though the oxygen abun- 
dance is less centrally peaked than for instan ce the iron abun- 
dance (IWerner et al.ll2006t ISimionescu et al ] 12009). There is a 
spread in abundance ratios across the different clusters of galax- 
ies studied, with for in stance the Ne/Fe ratio varying between 
0.73+0.18 in Hydra A (Simionescu et al. "2009) and 1.40±0.11 
in 2A 0335+096 (Werner etal. 2006). The best fit Fe/O ratio 
varies between 0.63 to 2.5 in a sample of six clusters of galaxies 
(Hydra A, Fornax, M 87, Centaurus, Sersic 159-03 and A 1060), 
and depends on the distance from the cluster core at which the 
measurement was made. 

The gas that enriches the ISM, i.e. that from stellar winds and 
SN explosions, is also expected to enrich the ICM. Therefore, 
naively one would expect that the metal abundance ratios mea- 
sured in the outflows of Seyfert 1 in galaxies and the ICM might 
be similar. The absolute abundances will differ because of the 
large reservoir of primordial gas in the ICM with which the en- 
riched gas will mix. A possible difference is expected if the com- 
bined star-formation history and initial mass function (IMF) of 
the galaxies belonging to the central part of the cluster is differ- 
ent from the star-formation history and IMF of Seyfert galax- 
ies. This would result in a different ratio of type la to core- 
collapse SNe, as well as alter the importance of stellar winds in 
the enrichment of the gas. The disturbed morphology of Mrk 509 
likely implies that it has had a different star-formation history 
from that of the elliptical galaxies dominating in the centres of 
clusters. 

In Table[3] we compare the O/Fe, Ne/Fe, Mg/Fe, Si/Fe, S/Fe, 
and Ca/Fe ratios for Mrk 509 and the cores of th e three clus- 
ters ofgakxiesFK^dra A (ISimionescu ( 



I et alj 

03 (Ide Plaa et alJl2006l) . and 2A 0335+096 (IWerner e"tani2006l) . 
The data we use are for Hydra A the results obtained from a 
fit to the inner 3', for Sersic 159-03 a fit to the inner 4', and 
for 2A 0335+096 also the inner 4' for the O/Fe and Ne/Fe, but 
the inner 3' for Si/Fe, S/Fe, and Ca/Fe. The oxygen and neon 
abundances are taken from the RGS fit, in the case of Hydra A 
the fit with ^ dem+lT, and gdem for the EPIC data are used 
(ISimionescu et al. 2009); for Sersic 159-03 and 2A 0335+096 
the results, are fr om the wdem model (Ide Plaa et al ]|200S 
IWemer et ai]|2006l) . For Mrk 509, we use the RGS abundances 
measured using model 3 and assuming that all velocity compo- 
nents have the same abundances (see Table [Tj), with the excep- 
tion of Si taken from the xabs modelling of the LETGS spec- 
trum. In clusters of galaxies, owing to the high temperature of 
the plasma and the extent of the source, broadening the emis- 
sion lines detected by the RGS with the brightness profile of the 
source, the abundances are generally measured using the EPIC 
instruments. Hence, the relative abundances are measured with 
respect to iron, which is well constrained with the EPIC instru- 
ments due to the presence of the Fe-L and Fe-K complexes. The 
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temperature is too high, even in cool-core clusters, to derive a 
reliable abundance for either nitrogen or carbon as these ele- 
ments are fully ionised. The oxygen and neon abundances we 
list were determined from RGS spectra, as neon is blended with 
the Fe-L complex, and oxygen is blended with a local oxygen 
absorption in the EPIC data. In the comparison below we use 
the abundance ratios determined for the cluster core, as the oxy- 
gen abundance is generally poorly constrained at outer annuli, 
where there is no good RGS data available. In general the O/Fe 
ratio increases with distance from the core, but for the other el- 
ements the abundance ratios are consistent with no change with 
distance. However, the absolute abundance does decrease with 
distance from the cluster core. 

The O/Fe, Ne/Fe, and Si/Fe abundance ratios are larger in 
Mrk 509 than for the cores of clusters, while the Mg/Fe, S/Fe 
and Ca/Fe are within the range measured in the cluster cores. All 
abundance ratios are consistent with at least the abundance ratio 
measured in one of the three clusters Hydra A, 2A 0335+096 
and Sersic 159-03. The higher oxygen, neon, and silicum abun- 
dances in the outflow of Mrk 509 reflects the greater importance 
to the enrichment of the ISM by both/either core-collapse SNe 
and/or stellar winds. The disturbed morphology of Mrk 509 im- 
plies that a higher core-collapse SN versus type la SN rate in 
the recent history of this galaxy is a likely explanation of the 
larger abundance ratios of O and Ne relative to Fe. Furthermore, 
stellar winds might be less likely to escape the potential of the 
galaxy, and therefore preferentially enrich the ISM, which would 
explain the overall higher abundance ratios observed in the out- 
flow in Mrk 509 than in the cores of clusters of galaxies. 



Table 3. Comparison of the relative abundances in Mrk 509 (see 
Table[T]i with those derived for the cores of the cluster of galaxies 
Hydra A, Sersic 159-03 and 2A 0335+096 (ISimionescu et alJ 
l2009llde PlaaetaPIIOOa IWerner et al.ll2006l) . 





Mrk 509 


Hydra A 


Sersic 159-03 


2A 0335+096 


O/Fe 


1.17±0.06 


0.76+0.11 


0.87+0.10 


0.49+0.05 


Ne/Fe 


1.31±0.17 


0.84+0.20 


0.73+0.12 


0.85+0.08 


Mg/Fe 


0.80±0.16 




0.32+0.04 


1.07+0.08 


Si/Fe 


1.53±0.60 


0.65+0.05 


0.69+0.02 


1.35+0.03 


S/Fe 


0.67±0.15 


0.58+0.05 


0.52+0.02 


1.14+0.03 


Ca/Fe 


1.05±0.26 


1.28+0.16 


1.04+0.05 


1.58+0.06 



5. Conclusions 

We have studied the relative abundances determined from the 
X-ray spectrum of the outflow in Mrk 509, which has allowed 
an accurate and reliable determination of the abundance ratio 
for five elements C/O, N/0, Ne/0, Mg/0, S/O, and Fe/O abun- 
dances. For Si/O and Ca/O, less accurate abundances were ob- 
tained. We have found that the measured relative abundances of 
C, N, Ne, Mg, Si, Ca, and Fe relative to oxygen, are consistent 
with proto-solar ratios. Sulphur is slig htly un derabundant rela- 
tive to the proto-solar ratios given bv iLodders & Palmd (120091) . 
These abundances have been confirmed using absorption mea- 
sure distribution modelling, albeit with larger errors. 

We compared our relative abundance determination to the 
absolute abundances measured for Mrk 279, another Seyfert 1 
galaxy at a similar redshift. For Mrk 279, the absolute abun- 
dances were determined from the HST-STIS and FUSE spec- 
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tra, and are consistent with the relative abundances measured 
for Mrk 509. The abundances derived for the warm absorber 
in IRAS 13349+2438 and NGC 3783 are consistent with our 
determined abundances and proto-solar ratios. We also com- 
pared our method and results with the abundance determination 
from the broad emission lines in quasars, and found the latter 
method to be unreliable because of the effect of the resonance- 
scattered Lya photons at the wavelength of the Nv broad- 
emission line. We therefore compared our results with the abun- 
dances measured for the Galactic centre and found that within 
the error bars and uncertainties the results are consistent. In ad- 
dition using the a bsolute abundances determined for Mrk 279 
(lArav et al.ll2007h . we compared the abundances within the nu- 
cleus of these Seyfert galaxies with ISM metallicities predicted 
from the mass-metallicity relationship for star-forming galax- 
ies (Tremonti et al. 2004). For Mrk 279, the absolute O abun- 
dance derived from UV absorption lines is consistent with the 
ISM metallicity, but for both measurements the error bars and 
uncertainties are large. As an alternative we compared our abun- 
dance ratios to those measured for the hot ISM in local elliptical 
galaxies. The Mg and Si abundance relative to Fe in Mrk 509 are 
consistent with those in the hot halo ISM, although we measure 
a larger O/Fe ratio than that of the hot ISM. 

Finally we compared the abundance ratios for O/Fe, Ne/Fe, 
Si/Fe, S/Fe, and Ca/Fe to the ratios determined in the cores of 
three clusters of galaxies. We found that the outflow has a higher 
O/Fe and Ne/Fe abundance than the cluster cores, while the 
Si/Fe, S/Fe, and Ca/Fe are within the range of measured clus- 
ter abundance ratios. This indicates that core-collapse SNe and 
stellar winds have been more important to the enrichment of the 
interstellar medium of Mrk 509 than to the enrichment of the 
intra-cluster medium of these clusters. Likely explanations are 
that stellar winds have not escaped the potential of the galaxy as 
readily as SNe ejecta, and that owing to a disturbance of the host 
galaxy of Mrk 509 the recent star-formation history is likely dif- 
ferent from that of cD galaxies residing at the cores of clusters. 
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